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ABSTRACT 


Leading charm hadrons produced in hadron-nucleus interactions cannot be 
adequately described within the parton fusion model. Recent results on 
charm baryon production in interactions at 330 GeV with the WA89 
detector disagree with fusion predictions. Intrinsic heavy quark pairs in the 
T,~{dds) wavefunction provide a simple mechanism for producing fast charm 
hadrons. We calculate leading charm baryon production from 7r“ and p 
projectiles in a two component model combining parton fusion with intrinsic 
charm. Final state D~, and A+ da/dxp distributions and D~/D^, 

D~/D'^ and A+/A^ asymmetries are compared to WA89 data. Predictions 
are made for 650 GeV and n~A interactions in the SELEX detector at 
Fermilab and for 800 GeV pA interactions. 
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1 Introduction 


One of the most striking features of charm hadroproduction is the leading particle effect: 
the strong correlation between the quantum numbers of the projectile and the final- 
state charm hadron. For example, more D~ than D'^ are produced at large xp in 
n~A —> D^X interactions B I, I, i i- Such correlations are remarkable because they 
explicitly contradict the perturbative QCD factorization theorem which predicts that 
heavy quarks hadronize through jet fragmentation functions independent of the initial 
state. 


While leading charm effects are well established for D mesons, observations of charm 
baryons are more rare |^. Two experiments with Il~{dds) beams promise to clarify the 
situation in the baryon sector. The hyperon beam, with a strange valence quark, presents 
a unique opportunity to study the flavor dependence of leading charm hadroproduction 
since both charm and charm-strange baryons can be leading. The first, WA89, which 
directs a 330 GeV hyperon beam on carbon and copper targets, has reported the xp 
distributions of D~(dc), T,^{ddc), S+(msc), and X^{udc) P as well as the D~/D~^, 
D ~/and Ac/Ac production asymmetries p. The second, SELEX |^, has a large 
acceptance for forward charm production, enhancing the charm baryon yield. Their 650 
GeV beam, approximately half 7T~{ud) and promises to improve current samples 
from both beams by up to an order of magnitude. They also plan to study the A 
dependence of leading charm. 


In previous work 00. a QGD mechanism which produces leading charm at large 
xp was introduced. An important feature of the model is coalescence, the process 
through which a charm quark hadronizes by combining with quarks of similar rapidities, 
such as projectile spectator valence quarks. In a gauge theory the strongest attraction is 
expected to occur when the spectators and the produced quarks have equal velocities . 
Thus the coalescence probability should be largest at small relative rapidity and rather 
low transverse momentum where the invariant mass of the Qq system is minimized, 
enhancing the binding amplitude. 


This coalescence occurs in the initial state where the projectile wavefunctions of e.g. 
the 7r“, p and can fluctuate into Fock conhgurations containing a cc pair such as 
\udcc), \uudcc) or \ddscc) respectively. In these states, two or more gluons are attached 
to the charm quarks, reducing the amplitude by 0{af) relative to parton fusion |^. The 
longest-lived fluctuations in states with invariant mass M have a lifetime of 0{2PiaP,/M‘^) 
in the target rest frame where Piab is the projectile momentum. Since the comoving 
charm and valence quarks have the same rapidity in these states, the heavy quarks 
carry a large fraction of the projectile momentum and can thus readily combine to 
produce charm hadrons with large longitudinal momentum. Such a mechanism can then 
dominate the hadroproduction rate at large xp. This is the underlying assumption of 


the intrinsic charm model in which the wavefunction fluctuations are initially far off 
shell. However, they materialize as charm hadrons when light spectator quarks in the 
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projectile Fock state interacts with the target |15]. Since such interactions are strong, 


charm production will occur primarily on the front face of the nucleus in the case of a 
nuclear target. Thus the intrinsic charm mechanism has a stronger A dependence than 
charm production by leading-twist fusion. 


In this work, we concentrate on the charm hadrons studied by WA89 and SELEX 
in order to further examine the relationship between fragmentation and coalescence 
mechanisms. The calculations are made within a two-component model: leading-twist 
fusion and intrinsic charm 
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Leading particle correlations are also an integral part of the Monte Carlo program 
PYTHIA based on the Lund string fragmentation model. In this model it is assumed 
that the heavy quarks are produced in the initial state with relatively small longitudinal 
momentum fractions by the leading twist fusion processes. In order to produce a strong 
leading particle effect at large xf, the string has to accelerate the heavy quark as it 
fragments into the hnal-state heavy hadron. Such a mechanism demands that large 
changes of the heavy quark momentum take place in the hnal state. Other models 
of leading charm production by recombination in the hnal-state have been suggested 
||T^ , |20[| . However, in this work we will only compare our results with the commonly 
used PYTHIA Monte Carlo. 


In this paper, we hrst discuss the conventional mechanism for charm production at 
leading twist, parton fusion, and how the hyperon beam is taken into account in the 
model. Section 3 reviews the intrinsic charm model and describes the extension of the 
model used in this work. In section 4, we compare our results on interactions with 
the WA89 data and make predictions for SELEX with E~ and 7r~ beams as well as pA 
interactions. Finally, we summarize our results. 


2 Leading-Twist Charm Production 


In this section we briehy review the conventional leading-twist model for the production 
of charm hadrons in E~p, pp and Ti~p interactions. In leading-twist QCD, heavy quarks 
are produced by the fusion subprocesses gg —> QQ and qq —> QQ. The factorization 
theorem ^ predicts that fragmentation is independent of the quantum numbers of both 
the projectile and target. We will also show the corresponding distributions of charm 
hadrons predicted by the PYTHIA model 


m 


Our calculations are at lowest order in A constant factor iP ~ 2 — 3 is included 
in the fusion cross section since the next-to-leading order xp distribution is larger than 
the leading order distribution by an approximately constant factor |^]. An additional 


factor of two is included to obtain the single charm distribution, twice the cc cross section. 
Note that neither leading order production nor the next-to-leading order corrections can 
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produce flavor correlations |^] such as those observed in leading charm production. 
The charm hadron xp distribution, where xp = {2mp/y/s) sinh?/, has the factorized 


form 17 


da 

dxp 


= —XaXb 


HAB{Xa,Xb) — 


1 Dh/c{z^) 


^3 


dz3 dy2 dp. 


T ) 


( 1 ) 


where a and b are the initial partons, 1 and 2 are the produced charm quarks with 
rric = 1.5 GeV, and 3 and 4 are the final-state charm hadrons. The convolution of the 
subprocess cross sections for qq annihilation and gluon fusion with the parton densities 
is included in HAB{xa,Xb), 


ddr, 


HABiXa,Xb) = + f^{Xa)fg{Xb)]-^ + fg{Xa)fg{Xb) 


da, 


99 


di 


( 2 ) 


where A and B are the interacting hadrons and the scale dependence of the parton 
densities has been suppressed. The subprocess cross sections can be found in Ref. |2^ . 
Since we study cc production with several different projectiles, we specify the general qq 
convolution for three light flavors: 


ff{Xa)fg{Xb)\ = U^{Xa)u^{Xb) + U^{Xa)U^{Xb) 

<? 

+d^{Xa)d^ {Xb) + d^{Xa)d^{Xb) + S^{Xa)s^{Xb) + S^{Xa)s^{Xb) (3) 


Parton distributions of the hyperon are not available. However, using baryon number 
and momentum sum rules, a set of parton distributions for the can be inferred from 
the proton distributions: 



u^{x) dx 


1 


dl {x) dx 


/ d^ (x)dx 
Jo 

/ sf (x) dx 
Jo 


2 

1 . 


( 4 ) 

( 5 ) 


We also identify s^(x) = {x). Similar relations can be made for the sea quarks. The 

gluon distributions are thus assumed to be the same in the and the proton. Both the 
GRV LO 94 [Q and MRS parton distribution functions with u ^ d ^ s were 

used. Other, older distributions with a symmetric sea, u = d = s, produce identical 
results for and pp interactions. 


The fragmentation functions, Dh/c{z), describe the hadronization of the charm quark 
where 2 ; = Ph/Pc is the fraction of the charm quark momentum carried by the charm 
hadron, assumed to be collinear to the charm quark. According to the factorization 
theorem, fragmentation is independent of the initial state and thus cannot produce 
flavor correlations, precluding a leading charm effect. This uncorrelated fragmentation 
will be modeled by two extremes: a delta function, 6{z — 1), and the Peterson function 
2^, extracted from e+e“ data. The Peterson function, derived from a non-hadronic 
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initial state, predicts a softer xp distribution than observed in hadroproduction, even at 
moderate xp since the fragmentation decelerates the charm quark. The parameters 
of the Peterson function we use here are taken from e+e” studies of D production . 
Typically hts to charm baryon fragmentation functions suggest increased deceleration of 
the charm quark in hnal-state baryons relative to mesons. On the other hand, the delta- 
function assumes that the charm quark coalesces with a low-x spectator sea quark or a 
low momentum secondary quark with little or no momentum loss [O. This assumption 

fl than Peterson 


is more consistent with low pp charm hadroproduction data pO 
fragmentation. 
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In Fig. |l|(a) we show the inclusive xp distributions calculated for both types of 
fragmentation in interactions at 330 GeV. Both sets of parton distributions are also 
shown. Very little difference in either total cross section or shape of the xp distributions 
can be discerned between the two sets of parton distributions. The delta function results 
in harder distributions than those predicted by Peterson fragmentation for xp > 0.2. 
However, as shown in [Q, even with this hard fragmentation the fusion model cannot 
account for the shape of leading charm baryon distributions. Figure ^(b) shows the 
relative contributions from gg fusion and qq annihilation to the total cross section at 
650 GeV, the energy of the SELEX beam, using the GRV LO 94 parton densities. Gluon 
fusion clearly dominates the production until xp ~ 0.6. We have checked that this is 
also true at the lower energy of the WA89 experiment, 330 GeV. 


We compare the distributions with those from pp and tt~p interactions at the 
same energies with our two choices of parton distributions in Figs. ^ and ^ respectively. 
Since the differences between the and pp xp distributions are rather small due to 
the dominance of gluon fusion, in Fig. ^(a) we show the ratio T,~p/pp at 330 GeV. The 
differences between and pp production are somewhat larger for the GRV LO 94 
distributions than the MRS D—', due to the relative assumptions of u/d. The GRV LO 
94 set is based on more recent data than the MRS D—' and should thus more accurately 
reflect the sea quark abundancies in the proton. In contrast, assumptions concerning 
charm quark fragmentation do not strongly affect the relative rates. In Fig. ^(b) we 
see that the relative qq contribution to pp production is somewhat larger at ~ 0 
than in production at the same energy but this does not affect the point where 
gg fusion ceases to dominate cc production. The pion valence distributions are harder, 
allowing charm production at larger xp than with a baryon beam, as shown in Fig. |^. 
For the pion, we use the GRV LO pion set with the GRV LO 94 proton set and with 
the MRS D—' distributions we use the SMRS P2 pion distributions |2^. However, the 
valence u quark in the ti~ does not change the relative importance of qq annihilation at 
650 GeV, as seen in Fig. ^(b). Much lower energies are needed for the vr” antiquark to 
lead to dominance of qq annihilation in cc production. 


The charm hadron distributions from PYTHIA 6.115 [jT^ for the three projectiles at 
330 GeV beam energy are also shown in Figs. |I|-^ (c) and (d). The PYTHIA calculations, 
based on 10^ events, use all default program settings along with the GRV LO 94 parton 
distributions. We note that in PYTHIA the hyperon valence quark distributions are 
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an average of the proton valence distribntions, = {u^ + d^)/3. In (c) the D~, 

Df, Ac and xp distribntions are shown while the D~, and distribntions 
are given in (d). The magnitnde of the cnrves reflect the relative abnndancies of charm 
hadrons prodnced by PYTHIA. 


The Lnnd string fragmentation model prodnces charm qnarks at string endpoints. 
The strings pnll the charm qnarks toward the opposite endpoints, typically the beam 
remnants. When the two string endpoints are moving in the same general direction, 
the charm hadron can be prodnced with larger longitndinal momentnm than the charm 
qnark. In the case where the string invariant mass is too small for mnltiple particle 
prodnction, a single hadron is prodnced [^, as in the S°(d(ic) and S°((isc) which share 
two valence qnarks with the These distribntions have a minimnm at xp ~ 0.3 
and 0.1 respectively and a peak at xp ~ 0.8, illnstrating the acceleration nndergone by 
charm quarks by strings with small invariant mass. The and are also accelerated 
by string fragmentation but the effect is not as strong because only one valence quark 
is in common with the projectile. 


In contrast, with a proton beam, as shown in Fig. g, only the Ac shows strong forward 
acceleration due to the common u and d quarks with the maximum in the xp distribution 
occurring at xp ~ 0.8. A second peak is notable for the but the acceleration effect 
is weaker for charm-strange baryon production, presumably due to the additional mass 
of the strange quark. While meson production does not show any signihcant leading 
behavior with a baryon projectile, the situation is reversed with the tt” beam where D 
and Dg production is clearly forward of all charm baryon production, produced centrally 
in the forward xp region. We will make further comparisons with PYTHIA when our 
full model is discussed. 


3 Intrinsic Particle Production 


The wavefunction of a hadron in QCD can be represented as a superposition of Fock 
state fluctuations, e.g. \nv), \nvg), \nvQQ), ...components where ny = dds for a 
uud for a proton and ud for a tt”. When the projectile scatters in the target, the 
coherence of the Fock components is broken and the fluctuations can hadronize either by 
uncorrelated fragmentation as for leading twist production or coalescence with spectator 
quarks in the wavefunction [M, The intrinsic heavy quark Fock components are 
generated by virtual interactions such as gg —>• QQ where the gluons couple to two or 
more projectile valence quarks. The probability to produce QQ fluctuations scales as 
oI^Mqq)/ rriQ relative to leading-twist production [Q. Intrinsic QQ Fock states are 
dominated by conhgurations with equal rapidity constituents so that, unlike sea quarks 
generated from a single parton, the intrinsic heavy quarks carry a large fraction of the 
parent momentum . 
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The frame-independent probability distribution of an n-particle cc Fock state is 


dP-r. 


dXi ■ ■ ■ dXr 




ml-Y.l=i{rnf/xi)y 


( 6 ) 


where iV„ normalizes the |ncc) probability, and n = 4, 5 for meson and baryon 
production from the \nvcc) conhguration. The delta function conserves longitudinal 
momentum. The dominant Fock configurations are closest to the light-cone energy shell 
and therefore the invariant mass, is minimized where mf = k‘^i + 

is the effective transverse mass of the particle and Xi is the light-cone momentum 
fraction. Assuming (k^j) is proportional to the square of the constituent quark mass, 
we choose niq = 0.45 GeV, = 0.71 GeV, and nic = 1-8 GeV 
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The intrinsic charm production cross section for a single charm hadron from the 
u-particle state can be related to P^ and the inelastic hN cross section by 


al(hN) = 


4m? 


(7) 


The factor of p^/4m^ arises from the soft interaction which breaks the coherence of the 
Fock state. To set the scale of the coherence factor /i we asssume that the NAS diffractive 
J/'ip cross section can be attributed to intrinsic charm. In this experiment the 
nuclear dependence oi J/'ip production in 7r“A interactions into a “hard” contribution 
with a nearly linear A dependence at low xp and a high xp “diffractive” contribution 
scaling as where /d = 0.77 for pion and 0.71 for proton beams, characteristic of soft 
interactions. Then we assume that the diffractive fraction of the production cross section 
331 is the same for charmonium and charm hadrons. In Ref. iTIIl, ~ 0.2 GeV^ was 


found, however, calculations with more recent parton densities suggest that jjF' ~ 0.1 
GeV^. We thus obtain (Tf^{TiN) k. 0.5 /rb and af^{pN) ^ 0.7 pb at 200 GeV. We take 
ps = 0.31%, as determined from an analysis of the EMG charm structure function data 
|3^ . A recent reanalysis of the EMG data with next-to-leading order calculations of 


leading twist and intrinsic charm electroproduction is consistent with the presence of an 
intrinsic charm component in the proton at large x-q-^ of ~ 1% or less |^. For simplicity. 


we will always assume that the total probability for a charm quark in an \nycc) state is 
0.31% pOl regardless of the projectile identity. 


The inelastic cross section has not been measured. However the total and elastic 
Ap cross sections have been parameterized for beam momenta less than 200 GeV albeit 
with large statistical uncertainties. Extrapolating these cross sections to 330 GeV, we 
found that <7“^ > a™ at this energy which seems unlikely. To be conservative, we 
therefore scaled a™ to at the highest measured Ap energy and used the energy 
dependence of cr™ thereafter to set the scale for E“p interactions at larger values of a/s. 

There are two ways of producing charm hadrons from intrinsic cc states. The first 
is by uncorrelated fragmentation, previously discussed in Section 2. Additionally, if the 
projectile has the corresponding valence quarks, the charm quark can also hadronize 
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by coalescence with the valence spectators. The coalescence mechanism thus introduces 
flavor correlations between the projectile and the flnal-state hadrons, producing e.g. S°’s 
with a large fraction of the momentum. 

First we briefly discuss charm production by uncorrelated fragmentation. If we as¬ 
sume that the c quark fragments into a D meson, the D distribution is 

J pnF . n j pn 

= / dzWdXi- -— Dd/c{z)5{xd - zx^) , ( 8 ) 

dXj;) J dXi . . . dXn 

These distributions are assumed for all intrinsic charm production by uncorrelated frag¬ 
mentation with Dh/c{z) = S{z — 1). We will not use Peterson function fragmentation 
further in this work. 


The coalescence distributions, on the other hand, are specific for the individual charm 
hadrons. It is reasonable to assume that the intrinsic charm Fock states are fragile and 
can easily materialize into charm hadrons in high-energy, low momentum transfer reac¬ 
tions through coalesence. These contributions, taken from Ref. []^, do not include any 
binding energy of the produced hadrons or any mass effects. The coalescence contribu¬ 
tion to charm hadron production is 

^ pnC p n j pn 

= j \{dxi— - ^^j—5{xh-xhi -• (9) 


dx 


H 


i=l 


' dxi . . . dXr, 


The coalescence function is simply a delta function combining the momentum fractions 
of the quarks in the Fock state configuration that make up the valence quarks of the 
flnal-state hadron. 


We now compare and contrast D~, , D~ and meson and and 

baryon production by coalescence from p and vr” projectiles. We note that 
not all of these hadrons can be produced from the minimal intrinsic charm Fock state 
configuration, |nycc). However, coalescence can also occur within higher fluctuations 
of the intrinsic charm Fock state. For example, in the proton, the and can be 
produced by coalescence from \nvccdd) and \nvccss) configurations. These higher Fock 
state probabilities can be obtained using earlier results []^ ^ . In a previous study of i/’V' 
production from \nvcccc) states [Q the double intrinsic charm production probability, 

arise from these 
tt^N) ^ 20 pb set 
This value can then be used to 


Pice, was determined assuming that all the measured pairs 
configurations. The resulting upper bound on the model, = Uj 
21 


by experiment 


requires Py 


4.4% Pc PI, m 


estimate the probability of light quark pairs in an intrinsic charm state. We expect that 
the probability of additional light quark pairs in the Fock states to be larger than pec. 


Pcq~ ^ Pcc, (10) 

X'FlqJ 

leading to Peu = Ped ~ 70.4% Py and Pcs ~ 28.5% Py. To go to still higher config¬ 
urations, e.g. for production from a 7i~, one can make the similar assumption that 
Pcsu = 70.4%Pcs. 


























Particle 

S (dds) 

piuud) 

71 iud) 

D-{dc) 

5 

5 

4 

D+{dc) 

7{dd) 

7{dd) 

6{dd) 

Ac{udc) 

7{uu) 

5 

6{uu) 


5 

7{dd) 

Q{dd) 

D;{sc) 

5 

7(ss) 

6(ss) 

Dtiac) 

7(ss) 

7(ss) 

6(ss) 

S^(dsc) 

5 

7(ss) 

6(ss) 

S+(msc) 

7{uu) 

7{ss) 

8{ssuu) 


Table 1: The lowest number of partons needed in an intrinsic charm Fock state con- 
hguration for the charm particle to be produced by coalescence. Note that 4 and 5 
correspond to the minimal |nycc) conhguration while the higher states refer to \nvccdd) 
etc. 


In Table we show the minimum number of partons needed in each conhguration 
to produce a given charm hadron. When more than the minimal |nycc) state is nec¬ 
essary for coalescence to occur, the additional light quark pairs required in the state 
are indicated. While we include the eight particle conhguration necessary to produce a 
by coalescence from a vr” projectile, we will conhne our discussion to charm hadron 
production from the minimal state and states with one additional qq pair only. 


The total intrinsic charm contribution to charm hadron production is a combination 
of uncorrelated fragmentation and coalescence. In previous works 0,PI |I6|, [TTH only 
production by uncorrelated fragmentation from the minimal |nycc) states and coales¬ 
cence from the minimum Fock state conhguration was considered. This was because a 
signihcant leading ehect is present only in the minimal conhguration, i.e. there is no 
diherence between and D~ mesons produced from \nv(xdd) states. Also, as more 
partons are included in the Fock state, the coalescence distributions soften and approach 
the fragmentation distributions, eventually producing charm hadrons with less momen¬ 
tum than uncorrelated fragmentation from the minimal cc state if a sufficient number of 
qq pairs are included. There is then no longer any advantage to introducing more light 
quark pairs into the conhguration—the relative probability will decrease while the po¬ 
tential gain in momentum is not signihcant. However, if some fraction of the hnal-state 
hadrons are assumed to be produced from higher Fock conhgurations, then all possible 
hnal-states from those conhgurations should also be included. Therefore in this paper, 
we consider production by fragmentation and coalescence from the minimal state and 
the next higher states with uu, dd and ss pairs. 


The probability distributions, {1/P^){dP^/dxu), are given in Figs. for p 
and 7r“ projectiles respectively. It is clear from Fig. ^ that the projectile allows the 
greatest coalescence production of charm hadrons from the minimal Fock conhguration. 
Fig. ^(a). The charm baryons are quite fast, taking more than 50% of the projectile 
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momentum. The difference between charm and charm-strange hadron production is 
very small due to the strange and light quark mass difference. Because the strange 
quark is more massive, it carries a somewhat larger fraction of the momentum than 
the light quarks, resulting in a slightly larger average momentum for the and the 
D~ relative to the and D~, on the order of 3-4% as can be seen in Table ^ The 
c quark distribution itself, leading to uncorrelated fragmentation, carries ~ 25% of the 
projectile momentum in the minimal state. This is reduced by 35% in the seven- 
particle Fock conhgurations. In this model, the c and c probability distributions are 
identical. We note that these higher conhgurations can produce, for example, and 
Ac baryons from \nvccdd) and \nvccuu) states respectively with the same probability 
distribution, shown in Fig. ib), but not necessarily with the same relative probability, 
as we will show shortly. Introducing an ss pair to the 7-particle conhguration reduces 
the average momentum of the hnal state hadron by ~ 2% over the average in the 7- 
particle conhgurations with lighter qq pairs. In addition to the reduction of the average 
momentum of the c quark in the higher conhgurations, the hnal-state charm hadron 
momentum from this conhguration is reduced as well, suggesting that no more signihcant 
contribution to the overall momentum of the hnal hadron will be obtained by including 
yet higher Fock conhgurations. 

While fewer charm hadrons can be directly produced from the minimal conhguration 
of a proton projectile, as evident from Fig. their average momentum is somewhat 
higher than the due to the abscence of the strange valence quark. However, this 
only ahects the hnal-state average momentum by 1-2%. Final-state charm hadrons from 
a pion projectile, shown in Fig. |] have, on average, 20% more momentum than from a 
baryon projectile because the total velocity is shared between fewer initial partons. Note 
also that mesons from a four-particle Fock conhguration and baryons from a six-particle 
Fock state each receive half of the projectile momentum. 


4 Model Predictions 


We now turn to specihc predictions of our model. We begin with the xp distribution of 
the hnal-state charm hadrons. The xp distribution for hnal-state hadron H is the sum 
of the leading-twist fusion and intrinsic charm components, 

daf^N _ , da^ 

dxp dxp dxp ’ 


where da^/dxp is related to dP^/dxp by 

^ ^ in dPu 
dxp dxp 


( 12 ) 


The probability distribution is the sum of all contributions from the \nvcc) and the 
\nvccqq) conhgurations with q = u, d, and s and includes uncorrelated fragmentation 
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State 

Particle 

S (ny = dds) 

pijiv = uud) 

TT (ny = ud) 

nycc) 

c 

0.251 

0.256 

0.308 

)■) 

D-{dc) 

0.41 

0.419 

0.499 

)■) 

Kciudc) 

- 

0.58 

- 

)■) 

TP^ddc) 

0.573 

- 

- 

)■) 

B-(sc) 

0.427 

- 

- 

)■) 


0.59 

- 

- 

\nvccqq) 

c 

0.185 

0.188 

0.219 

)■) 

D-{dc) = D^{dc) 

0.31 

0.314 

0.359 

5!) 

Kciudc) 

0.433 

0.438 

0.5 

55 

TP^ddc) 

0.433 

0.438 

0.5 

55 

Djisc) 

0.32 

- 

- 

55 


0.444 

- 

- 

\nvccss) 

c 

0.179 

0.181 

0.211 

55 

D-{dc) 

0.302 

0.306 

0.349 

55 

Kciudc) 

- 

0.429 

- 

55 

TP^ddc) 

0.424 

- 

- 

55 

D-{sc) = Dt{sc) 

0.312 

0.316 

0.361 

55 

S^(dsc) 

0.434 

0.439 

0.5 

55 

S+(msc) 

0.434 

- 

- 


Table 2: The average value oi xp for charm particles produced by coalescence from S , 
p and 7r“ projectiles in \nvcc), \nyccqq) and Inyccss) states. In this case, qq = uu, dd. 
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and coalescence when appropriate, as described below. We use the same fragmentation 
function, either the delta or Peterson function, to calculate uncorrelated fragmentation 
in both leading twist fusion and intrinsic charm. In this section, we use only the delta 
function. 

Since experimental information on the relative rates of charm hadron production 
is incomplete, we assume that all the lowest lying charm hadrons produced by uncor¬ 
related fragmentation have equal probability in both leading-twist fusion and intrinsic 
charm. There are 10 charm hadrons—and the same number of anticharm hadrons—if 
excited charm hadrons such as D* and A+(2593) are excluded. Therefore the probabil¬ 
ity distribution for uncorrelated fragmentation into each of these hadrons is 10% of the 
total probability. As can be seen in Table 1, only a fraction of the possible hnal-state 
hadrons can be produced by coalescence. We use a simple counting scheme to arrive at 
the coalescence probability which enhances the production of leading charm at large xp- 
We note that the combined probability of fragmentation and coalescence of all charm 
hadrons cannot exceed the total production probability of the Fock state conhguration. 
Thus when a particular hnal-state hadron can be produced both by uncorrelated frag¬ 
mentation and coalescence, we multiply the sum of the fragmentation and coalescence 
probabilities by 0.5 to keep the total probability hxed. 

As a concrete example of how the total probability distributions of charm hadron 
production from the intrinsic charm model is calculated, we will describe and D~ 
production from the beam in our model in detail. The full complement of equations 
for all the hnal-state charm hadrons from p and vr” projectiles considered in this 
work can be found in the appendix. In the \ddscc) conhguration, there are four hnal- 
state hadrons with a valence c quark (2S°, and J/i/j) and also four hnal-state hadrons 
with a valence c quark {2D~, D~ and J/'ip)- Note that the J/'0 has been counted in each 
category. The D~ is then produced by coalescence with 50% of the total coalescence 
probability for hadrons with a valence c as well as by uncorrelated fragmentation of the 
c while the is only produced by uncorrelated fragmentation from this state. The 
probability distributions from this minimal Fock conhguration are then 

dPh- ^ 1 ( 1 ^ IdPifA 

dxp 2 yiO dxp 2 dxp ) 

dPl* _ 1 dPt” 

dxp 10 dxp 

where F refers to uncorrelated fragmentation and C to coalescence into the specihc hnal- 
state with the associated probability distribution, shown in Fig. H(a). The \ddsccqq) 
conhgurations where qq = uu, dd and ss allow coalescence production of eight hnal- 
state hadrons with a valence c and hve hnal-state hadrons with a valence c in each 
case. We will discuss and D~ from each of these conhgurations in turn. When 
qq = uu, the possible hadrons produced by coalescence are: 2S+, 2A+, and 

J/'ip with a valence c and 2D~, D~, and J/'0 with a valence c. A hnal-state D~^ 
can be produced by coalescence from the dd conhguration in one of the eight possible 


(13) 

(14) 
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final-state hadrons with a valence c quark (3S°, 3S°, and J/^j) while the D~ is 
produced by coalescence in three out of hve combinations {3D~, D~ and J/ip) with a 
valence c. The ss conhguration yields no by coalescence—4S°, r2°(ssc), Df 

and J/ip are allowed—while 2D~ are allowed out of 5 possible hadrons with valence c 
quarks— 2D~, 2D~ and J/i^. Finally, the total intrinsic charm probability distribution 
for these mesons is: 



dP£)+ 

dxp 


1 dP!f 1 

+ I + - 


1 dPE ^ 2dPl^ 


2 yiO dxp 2 dxp ) 2 yiO dxp 5 dxp ) 


1 dPS 


+ 


2 yiO dxp 5 dx 

L‘1PL + PPK+ 1 

10 dxp 10 dxp 2 


3dPS\ 1 2dPlS\ 


2 \^10 dxp 5 dxp J 


1 dPilf ^ ldP^,^\ 


+ 


1 dPif 


10 dxp 8 dxp / 10 dxp 


(15) 

(16) 


Lastly we note that only fragmentation from the minimal Fock state was included along 
with coalescence from the lowest possible conhguration in earlier work |^, |^. This 

corresponds to the hrst term of the D~ probability distribution while the D~^ distribution 
would be proportional to 0.5((l/10)dPj^^/(ia;i? -|- {l/8)dP^^^/dxp). 


We must also account for the fact that most of the data are taken on nuclear targets. 
In this case, the model assumes a linear A dependence for leading-twist fusion and an 
A°‘ dependence for the intrinsic charm component ||^ 

0.71 for protons (and S“) 


where a = 0.77 for pions and 


dc^hA 

dxp 


A d^\t 
A-AL 

dx jp 


+ 71" 


d^ 

dxp 


(17) 


This A dependence is included in the calculations. The intrinsic charm contribution to 
the longitudinal momentum distributions per nucleon is thus reduced for nuclear targets. 


We now compare the model calculations, both the full model of eqs. (|^ and 
and the simpler model used previously, to the WA89 data [|] on carbon and copper 
targets in Fig. ^ Since the data are unnormalized, we have normalized all curves to 
the first data point. The dot-dashed and dotted curves are results with the previous 
simplihed model |n, 00 on carbon and copper targets respectively. The full model is 
illustrated in the solid and dashed curves for the same targets. The agreement with the 
data is quite reasonable given both the low statistics of the data and our normalization 
to the hrst data point rather than htting the normalization to the data. The differences 
in the model distributions are most obvious for the shown in Fig. 0(c). The simpler 
model emphasizes the coalescence production from the \nvcc) state only. As can be seen 
from Fig. 0 and Table 0, the average of the coalescence distribution is more than a 
factor of two larger than that of a production by independent fragmentation of a c 
quark, producing a shoulder in the xp distribution, particularly for the carbon target 
(dot-dashed curve). Because the is produced by coalescence with 30% less average 
momentum from the 7-particle Fock states, the intermediate xp region is partially hlled 
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This Model 

Particle 

E-Gu (330) 

E-Gu (650) 

TT Gu (650) 

pGu (800) 

D-{dc) 

0.192 

0.147 

0.169 

0.120 

D+ldc) 

0.133 

0.112 

0.152 

0.107 

Ac{udc) 

0.145 

0.118 

0.154 

0.146 


0.187 

0.140 

0.154 

0.107 

D;(sc) 

0.165 

0.129 

0.151 

0.106 

DtM 

0.132 

0.111 

0.151 

0.105 


0.221 

0.160 

0.151 

0.106 

S+(nsc) 

0.160 

0.126 

0.150 

0.106 

PYTHIA 

Particle 

S-p (330) 

E-p (650) 

TT p (650) 

pp (800) 

D-{dc) 

0.14 

0.126 

0.254 

0.113 

D+{dc) 

0.18 

0.159 

0.173 

0.160 

Aciudc) 

0.54 

0.468 

0.153 

0.604 

TP^ddc) 

0.72 

0.707 

0.146 

0.35 

D-{sc) 

0.155 

0.139 

0.172 

0.097 

DtCsc) 

0.171 

0.154 

0.16 

0.153 


0.76 

0.767 

0.157 

0.123 

S+(nsc) 

0.55 

0.477 

0.156 

0.155 


Table 3: The average value of for charm particles produced in the full model for 
p and 7i~ projectiles on a copper target. The model results are given at 330 GeV and 
650 GeV for S“Gu interactions, 650 GeV for 7r“Gu interactions and 800 GeV for pGu 
interactions. For comparison, the average from PYTHIA at > 0 in each case is also 
shown. 

in, resulting in a smoother xp distribution even though the probability is reduced for 
the higher Fock states. Similar results can be seen for the other charm hadrons in Fig. |^. 

Note that the model results are in much better agreement with the data than the 
PYTHIA simulations at the same energy with the default settings, shown in Fig. |l](c) and 
(d). If the PYTHIA predictions are superimposed on Fig. ^ with the same normalization 
as our model, the PYTHIA results would considerably exceed the data at large xp for 
the charm baryons. In particular, since the c quark is pulled forward by a valence dd 
diquark, the rate from PYTHIA at xp ~ 0.8 would exceed the data by nearly four 
orders of magnitude. The differences in the results are also obvious in Table ^ where the 
average of all the model distributions on a copper target are compared to PYTHIA 
calculations with a proton target at the same energy. 

Another way to quantify leading charm production is through the asymmetry be- 
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tween leading and nonleading charm. The asymmetry is dehned as 


A{xf) = 


daildxF — daNh/dxF 
dai/dxF + dajsiL/dxF 


( 18 ) 


where L represents the leading and NL the nonleading charm hadron. High statistics 
data has previously been available only from vr” beams where a signihcant enhancement 
of D~ over production was seen at xf > 0.3 in qualitative agreement with 

]I^ also correctly predicted the 


111. The model 


the intrinsic charm calculation of Ref. 
symmetric production of D~ and mesons and A+ and 

E, SEl- 


baryons by vr beams 


Statistics are unfortunately limited on charm production by baryon beams. Recently 
the WA89 collaboration has presented the D~/D^, D ~/and A+/A)^ asymmetries 
from their data |^. In Fig. we compare our calculations with both models to 
this data as well as show a prediction for the asymmetry between the D~ and 
both of which are produced from the partons of the minimal Fock conhguration. The 
full model gives a larger asymmetry betwen D~ and at low xf than the simpler 
assumptions of previous work O, O because D~ production at intermediate xf is 


enhanced by coalescence production from the 7-particle conhgurations, see also Fig. |^(a). 
Our results with the full model are in qualitative agreement with the data, shown in 
Fig. ^(a). The measured D~/Df and A+/A)^ asymmetries are larger than our predictions 
at intermediate xf- The probability distribution for A)^ in our model is 


dPr 

dXF 


1 dPjf ^ 1 dPj! 

10 dxF 10 dxF 10 dxF 10 dxF 


(19) 


Some of the discrepancies between the model and the data may arise from the relatively 
low statistics of the Dg and Ac measurements. Our model is also quite crude in overall 
normalization for the different hnal states since we assume that all hnal-state hadrons 
are produced by independent fragmentation with the same probability. Not enough 
high statistics data exist yet for us to use experimental absolute production rates as 
a guide. The asymmetry between D~ and is interesting because the \ddscc) state 
of the can be thought of as a virtual D~E^ fluctuation, as has been suggested 
for proton fluctuations into K^A ^ and D~A'^ The is positive at 

hrst since the D~ xf distribution is larger at intermediate xf, especially when the 7- 
particle conhgurations are included. At larger xf, the baryon distributions always lead 
over charm mesons produced in the same conhguration, causing the asymmetry 

to approach —1 as xf —1- Similar results should be expected from the models of 
Refs. 


We now turn to predictions of charm hadron production at SELEX with 650 GeV 
beams of and tt" - m- First we give the charm hadron xf distributions for S Cu 
interactions and the relevant asymmetries in Fig. Since the leading-twist fusion cross 
section grows faster than the average xf of the particles studied decreases ~ 30% 
from 330 GeV to 650 GeV. A smaller decrease is found from the PYTHIA model, showing 
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the relative strength of the string fragmentation mechanism, as can be seen in Table |^. 
The is clearly the hardest distribution, followed by the The leads the 
because the more massive valence s quark carries more of the velocity than the d 
valence quarks. The leads the A+ in the 7-particle uu state for the same reason. The 
D~ and Dj, also produced from the 5-particle state have the hardest meson distributions 
but lag the baryons. The and have the softest distributions with the slightly 
harder because the quarks in the dd configuration get slightly more velocity than the 
ss conhguration with the more massive strange quarks. The asymmetries, which should 
be compared to the dashed curves in Fig. are somewhat reduced at higher energies, 
again due to the larger leading-twist cross section. 


Since SELEX will also measure charm hadroproduction with a tt beam at the same 
energy, these predictions are shown in Fig. 0 Because only the D is produced from 
the minimal Fock state conhguration, it shows the hardest xp distribution in Fig. pT|(a). 
Note that since cc production by leading-twist fusion alone is already signihcantly harder 
than the equivalent production by baryon projectiles, the distributions produced by 
coalescence from 6-particle conhgurations are not substantially enhanced over the fusion 
cross section, even at large xp- The intrinsic charm cross section is proportional to 
which increases slowly compared to the leading-twist cross section, further decreasing 
the predicted leading charm enhancement. Additionally, we note that charm baryons 
lead mesons produced by coalescence only in the 6-particle conhgurations since the 
baryons take ~ 50% of the pion momentum while the mesons take less, as seen in Fig. |] 
and Table |^. The PYTHIA distributions in Fig. aside from the leading D~, are 
more central, also evident from Table |^. In Fig. pT|(c), only the D~ /asymmetry is 
shown because the model predicts identical D~ and D'^ meson and Ac and A^ baryon 
distributions, see the appendix, hence no asymmetry. We note that the asymmetry is 
reduced compared to calculations at lower energy . 


The primary proton beam for hxed-target experiments at Fermilab is 800 GeV so 
for completeness, we also give predictions for a possible pA measurement at this energy 
in Fig. In this case, the A+ has the hardest xp distribution followed by the D~, 
both of which are produced by coalescence from the 5-particle Fock state. Again, the 
and are somewhat harder than the 5+ and distributions respectively due to 
the relative partitioning of the parton velocity in the 7-particle uu and dd conhgurations 
compared to the 7-particle ss state. The model predicts a strong D ~/asymmetry 
as well as a D~/ asymmetry, comparable to the D~ asymmetry predicted for 
the interactions. On the other hand, the D~/D^ asymmetry is quite weak. Such 
measurements with a proton beam would provide a useful complement to a high statistics 
measurement. A comprehensive understanding of data with proton projectiles has 
suffered in the past from a lack of statistics and high precision proton data, compared 
to that from and vr” projectiles, could eliminate certain models. 
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5 Summary and Conclusions 


We have refined the intrinsic charm model of Refs. 00.0. inclnding both the 
minimal Fock state and all the configurations with an additional qq pair. We have 
applied a simple counting scheme to determine the relative contribution of each state 
to the final charm hadron distribution. The model compares rather favorably to the 
xp distributions measured by WA89 p and produces reasonable agreement with their 
measured D ~/asymmetry while falling short of the D ~/and Ac/Ac data [Q at 
intermediate xp. 


Further, we have made predictions for charm hadron production at the energy of 
SELEX for both and ti~ projectiles. Predictions for production by an 800 GeV proton 
beam are also given. High statistics data on charm production from a combination 
of these projectiles could eliminate certain classes of models and perhaps distinguish 
between coalescence in the initial state, as in the intrinsic charm model, and in the final 
state, as in models such as PYTHIA []^|. The simple counting scheme employed here 


could be replaced with relative rates from data. However, the shapes of the distributions 
would not change significantly in our model. Therefore a collection of charm production 
data could define the role of intrinsic charm in future experiments. 


Acknowledgments: We thank S.J. 
Sjostrand for discussions. 


Brodsky, E. Ramberg, H.-W. Siebert and T. 


References 

[1] E.M. Aitala et al. (E791 Collab.), Phys. Lett. B371 (1996) 157. 

[2] M. Aguilar-Benitez et al. (LEBC-EHS Collab.), Phys. Lett. 161B (1985) 400; Z. 
Phys. C31 (1986) 491. 

[3] S. Barlag et al. (ACCMOR Collab.), Z. Phys. C49 (1991) 555. 

[4] M.I. Adamovich et al. (WA82 Collab.), Phys. Lett. B305 (1993) 402. 

[5] G.A. Alves et al. (E769 Collab.), Phys. Rev. Lett. 72 (1994) 812. 

[6] J.C. Collins, D.E. Soper and G. Sterman, Perturbative QCD, ed. A. H. Mueller 
(World Scientific, Singapore, 1989). G. Bodwin, Phys. Rev. D31 (1985) 2616, D34 
(1986) 3932. J. Qiu and G. Sterman, Nucl. Phys. B353 (1991) 105,137. 

[7] E.M. Aitala et al. (E791 Collab.), Phys. Lett. B411 (1997) 230. 

[8] R. Werding (WA89 Collab.), in Proceedings of ICHEP94, 27*^ International Con¬ 
ference on High Energy Physics, Glasgow, Scotland (1994). 


17 






[9] M.I. Adamovich et al. (WA89 Collab.), CERN-EP/98-41. 

[10] E. Ramberg in Hyperons, Charm and Beauty Hadrons, Proceedings of the 2'^'^ Inter¬ 
national Conference on Hyperons, Charm and Beauty Hadrons, Montreal, Canada, 
1996, edited by C.S. Kalman et al, Nucl. Phys. B (Proc. Suppl.) 55A (1997) 173. 

[11] R. Vogt and S.J. Brodsky, Nucl. Phys. B438 (1995) 261. 

[12] R. Vogt and S.J. Brodsky, Nucl. Phys. B478 (1996) 311. 

[13] S.J. Brodsky, J.F. Gunion and D.E. Soper, Phys. Rev. D36 (1987) 2710. 

[14] S.J. Brodsky, P. Hoyer, C. Peterson and N. Sakai, Phys. Lett. B93 (1980) 451; S. 
J. Brodsky, C. Peterson and N. Sakai, Phys. Rev. D23 (1981) 2745. 

[15] S.J. Brodsky, P. Hoyer, A.H. Mueller and W.-K. Tang, Nucl. Phys. B369 (1992) 
519. 

[16] R. Vogt, S.J. Brodsky and P. Hoyer, Nucl. Phys. B360 (1991) 67. 

[17] R. Vogt, S.J. Brodsky and P. Hoyer, Nucl. Phys. B383 (1992) 643. 

[18] T. Sjostrand, Comput. Phys. Commun. 82 (1994) 74. Program updates and docu¬ 
mentation can be found at \http://www. they, lu. se/tf2/statf'/torbjorn/Pythia. htm\ . 

[19] R.C. Hwa, Phys. Rev. D27 (1983) 653; Oregon preprint OITS-539 (1994). 

[20] V. A. Bednyakov, Dubna preprint JINR E2-94-79 (1994), |hep-ph/9403270| ; V.G. 
Kartvelishvili, A.K. Likhoded, and S.R. Slobospitskii, Sov. J. Nucl. Phys. 33(3) 
(1981) 434 [Yad. Fiz. 33 (1981) 832]. 

[21] R. Vogt, Z. Phys. C71 (1996) 475. 

[22] S. Frixione, M.L. Mangano, P. Nason, and G. Ridolh, Nucl. Phys. B431 (1994) 453. 

[23] R.K. Ellis, in Physics at the 100 GeV Scale, Proceedings of the 17*'^ SLAG Summer 
Institute, Stanford, California, 1989, edited by E.C. Brennan (SLAG Report No. 
361) 45. 

[24] M. Gliick, E. Reya, A. Vogt, Z. Phys. C67 (1995) 433. 

[25] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Lett. B306 (1993) 145. 

[26] M. Gliick, E. Reya, and A. Vogt, Z. Phys. C53 (1992) 651. 

[27] P.J. Sutton, A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D45 (1992) 
2349. 

[28] C. Peterson, D. Schlatter, 1. Schmitt, and P. Zerwas, Phys. Rev. D27 (1983) 105. 


18 






[29] J. Chirn, in proceedings of the ‘International Symposium on the Production and 
Decay of Heavy Flavors’, Stanford, CA, edited by E. Bloom and A. Fridman, (1987) 
131. 

[30] M. Aguilar-Benitez et al. (LEBC-EHS Collab.), Phys. Lett. B189 (1987) 476. 

[31] M. Aguilar-Benitez et al. (LEBC-EHS Collab.), Phys. Lett. B161 (1985) 400; Z. 
Phys. C31 (1986) 491. 

[32] S. Barlag et al. (ACCMOR Collab.), Phys. Lett. B247 (1990) 113. 

[33] T. Sjostrand, private communication. 

[34] P. Hoyer and S.J. Brodsky, in proceedings of the ‘Topical Conference on Particle 
Production Near Threshold’, Nashville, IN, edited by H. Hann and E.J. Stephenson 
(1990) 238. 

[35] J. Badier et al. (NA3 Collab.), Z. Phys. C20 (1983) 101. 

[36] J.J. Aubert et al. (EMC Collab.), Phys. Lett. IlOB (1982) 73; E. Hoffmann and 
R. Moore, Z. Phys. C20 (1983) 71. 

[37] B.W. Harris, J. Smith, and R. Vogt, Nucl. Phys. B461 (1996) 181. 

[38] R. Vogt and S.J. Brodsky, Phys. Lett. B349 (1995) 569. 

[39] J. Badier et al. (NA3 Collab.), Phys. Lett. 114B (1982) 457. 

[40] J. Badier et al. (NA3 Collab.), Phys. Lett. 158B (1985) 85. 

[41] R. Vogt, Nucl. Phys. B446, (1995) 159. 

[42] M.I. Adamovich et al. (BEATRICE Collab.), Nucl. Phys. B495 (1997) 3. 

[43] B.-Q. Ma and S.J. Brodsky, in proceedings of the ICTP Conference on Perspectives 
in Hadronic Physics, Trieste, 1997, |hep-ph/9707408 . 

[44] S. Paiva et al, IFUSP-P-1240, |hep-ph/9610310 . 


19 








Appendix 


Here we give the probability distributions for D~, D^, D~ and Df mesons and A^, 
and baryons for production by the minimal and hrst three higher Fock state 
conhgurations from proton and 7i~ projectiles. The probability distributions for 
each hnal state are given in Figs. |^-P. We note that the predictions for A+ and are 
identical in all cases because their quark content is the same. Recall that = 0.31%, 
PL = P-L = 70.4% pf and P-l, = 28.5% Pf. 


We begin with the In the \ddscc) conhguration, there are four hnal-state hadrons 
with a valence c quark (2S°, and J/t/’) and also four hnal-state hadrons with a valence 
c quark {2D~, D~ and J/V’). The \ddsccqq) conhgurations where qq = uu, dd and ss 
allow coalescence production of eight possible hnal-state hadrons with a valence c and 
hve possible hnal-state hadrons with a valence c. When qq = uu, the possible hadrons 
produced by coalescence are 2S+, 2A+, and with a valence c and 2D~, 

D~, and J/'ip with a valence c. The dd conhguration allows coalescence production 

of the following hadrons with a valence c quark, 3S°, 3S°, D~^ and J/ip, and, with a 
valence c, 3Z1~, Dp and J/ip- The ss conhguration yields 4S°, f2°(ssc). Dp and J/p) 

while the hnal-state valence c quarks hadrons are 2D~, 2Dp and J/p!. We have: 


dPo- 

dxp 


dPQ+ 

dxp 

dxp 

dxp 


dxp 


dPj^+ 

dxp 

dP^o 

_'~'c 

dxp 


1 1 nys; ‘idPL\ 

2 yiO dxp 2 dxp ) 2 yiO dxp 5 dxp j 

1 nys 3ys\ 1 /yyg 

2 yiO dxp 5 dxp j 2 yiO dxp 5 dxp j 

1 dPL . yyg , 1 ('lyg , lygj, 1 dPL 

10 dxp 10 dxp 2 yiO dxp 8 dxp j 10 dxp 

llA! + i ('llAI + llSSj + llAi + 1 dPL 

10 dxp 2 yiO dxp 4 dxp j 10 dxp 10 dxp 

1 IdPm 1 IdPgj 

2 yiO dxp 4 dxp j 2 yiO dxp 8 dxp j 

1 /yyg 3dPL\ 1 /yyg; i 

2 yiO dxp 8 dxp j 2 yiO dxp 8 dxp j 

1 /y yg iyg\ i /y yg i yg'i 

2 yiO dxp 4 dxp j 2 yiO dxp 5 dxp j 

1 /y yg iyg\ i g yg 2 yg\ 

2 yiO dxp 5 dxp j 2 yiO dxp 5 dxp j 

y yg y yg y yg i g yg lygA 

10 dxp 10 dxp 10 dxp 2 ylO dxp 8 dxp j 

1 /yyg lygA i /yyg i yg'i 

2 yiO dxp 2 dxp j 2 yiO dxp 8 dxp j 


(A-1) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 

(A-6) 
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(A-7) 


dP^+ 

_ 

dxp 


1 1 n ^ 1 

2 yiO dxp 8 dxp ) 2 yiO dxp 2 dxp ) 

1 iP? 1 }_dPS\ 1 

10 dxp 2 yiO dxp 4 dxp j 10 dxp 10 dxp 


(A-8) 


Fewer charm hadrons are produced by coalescence from the hve-quark conhguration 
of the proton since it has no valence strange quark. In the \uudcc) conhguration, there 
are four hnal-state hadrons with a valence c quark (2A+, S++(mmc) and J/'ip) and also 
four hnal-state hadrons with a valence c quark {2D , D~ and J/ip). The \uudccqq) 
conhgurations allow coalescence production of eight possible hnal-state hadrons with a 
valence c and hve possible hnal-state hadrons with a valence c. When qq = uu, the 
possible hadrons produced by coalescence are: 3A+, 3S++, and J /ip with a valence 
c and D~, 3D^ and J/ip with a valence c. The dd conhguration allows coalescence 
production of the following hadrons with a valence c quark, 4A+, S°, D~^ and J/t/’, 

and, with a valence c, 2D~, 21^ and J/ip. The ss conhguration yields 2S+, 2A+, 

S++, Dp and J /ip while the hnal-state valence c quarks hadrons are 2D~ ^ D~ and 
J/ip. Then: 


dPp,- 

dxp 


dPp+ 
dxp 

dxp 


dxp 

dPo- 

dxp 

dPj^+ 

dxp 

dP’^O 

_ "~~'c 

dxp 

dP^+ 

_ 

dxp 


1 ^ 1 dP«^\ 1 1 

2 yiO dxp 4 dxp j 2 yiO dxp 5 dxp j 

1 n ^ ■idPS\ 1 /"i^ 

2 yiO dxp 5 dxp j 2 yiO dxp 5 dxp j 

1 dP^ J_d^ 1 (]_d^ l^'l 

10 dxp 10 dxp 2 yiO dxp 8 dxp j 10 dxp 

1 1 dPS^\ 1 / 

2 yiO dxp 2 dxp j 2 yiO dxp 8 dxp j 

+1 (LiPl + + i ILPS. + D2SP 

2 yiO dxp 2 dxp j 2 yiO dxp 4 dxp j 
1 dpff 1 (]_^ 1 

10 dxp 10 dxp 2 yiO dxp 8 dxp j 10 dxp 


1 dP,f 1 

10 dxp 10 dxp 10 dxp 2 
1 df-^/ 1 
10 dxp 10 dxp 10 dxp 2 

1 dPff ^ 1 

10 dxp 10 dxp 10 dxp 2 
1 dP^ ^ ^ 1 
10 dxp 10 dxp 10 dxp 2 


1 dP"- ^ 1 dpm 

10 dxp 5 dxp J 


1 dP^ ^ 1 dPS\ 

10 dxp 8 dxp J 


1 dPi[ ^ 1 

10 dxp 8 dxp J 


1 dP^ ^ 1 dPS\ 

10 dxp 4 dxp ) 


(A-9) 

(A-10) 

(A-11) 

(A-12) 

(A-13) 

(A-14) 

(A-15) 

(A-16) 


Charm and anticharm hadron production is more symmetric from the tt because the 
projectile contains a valence antiquark of its own. In the minimal \udcc) conhguration. 
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there are two possible final-state hadrons with a valence c quark {D^ and J/t/’) and also 
two possible final-state hadrons with a valence c quark {D~ and J/'ip)- The \udccqq) 
configurations allow coalescence production of four possible final-state hadrons with a 
valence c and likewise four possible final-state hadrons with a valence c. When qq = uu, 
the possible hadrons produced by coalescence are: A+, 2D^ and J/'ip with a valence c 
and D~, and J/'ip with a valence c. The dd configuration allows coalescence 

production of the following hadrons with a valence c quark, and J/t/’, and, 

with a valence c, D~, and J/t/’- The ss configuration yields Df, and J/'ip 

with a valence c while the possible final-state valence c quarks hadrons are D~^ Dj, 
and J/ip. In this case: 
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Figure 1: Charm production by leading-twist fusion in interactions, (a) Two 

parton distribution functions with two different fragmentation functions are shown at 
330 GeV. The curves show calculations with the MRS D—' parton distributions with 
delta function fragmentation (solid) and the Peterson function (dot-dashed) and with 
the GRV LO 94 parton distributions with delta function fragmentation (dashed) and the 
Peterson function (dotted). In (b) calculations with the GRV LO 94 parton distributions 
with delta function fragmentation are given at 650 GeV for the qq component (dashed), 
gg component (dot-dashed) and the total production cross section (solid). Charm hadron 
production in PYTHIA 6.115 at 330 GeV is shown in (c) and (d) with the distributions 
labeled as indicated. 
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Figure 2: Charm production by leading-twist fusion in pp interactions, (a) The cross 
section ratios app/<JY,-p are given for two parton distribution functions with two different 
fragmentation functions at 330 GeV. The curves show calculations with the MRS D—' 
parton distributions with delta function fragmentation (solid) and the Peterson function 
(dot-dashed) and with the GRV LO 94 parton distributions with delta function fragmen¬ 
tation (dashed) and the Peterson function (dotted). In (b) calculations with the GRV 
LO 94 parton distributions with delta function fragmentation are given at 650 GeV for 
the qq component (dashed), gg component (dot-dashed) and the total production cross 
section (solid). Charm hadron production in PYTHIA 6.115 at 330 GeV is shown in (c) 
and (d) with the distributions labeled as indicated. 
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Figure 3: Charm production by leading-twist fusion in 'K~p interactions, (a) Two 
parton distribution functions with two different fragmentation functions are shown at 
330 GeV. The curves show calculations with the MRS D—' parton distributions with 
delta function fragmentation (solid) and the Peterson function (dot-dashed) and with 
the GRV LO 94 parton distributions with delta function fragmentation (dashed) and the 
Peterson function (dotted). In (b) calculations with the GRV LO 94 parton distributions 
using delta function fragmentation are given at 650 GeV for the qq component (dashed), 
gg component (dot-dashed) and the total production cross section (solid). Charm hadron 
production in PYTHIA 6.115 at 330 GeV is shown in (c) and (d) with the distributions 
labeled as indicated. 
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Figure 4: Charm hadron production in the intrinsic charm model with a projectile. 
The probability distributions, {1/P^){dP^/dxH), for uncorrelated fragmentation and 
coalescence are given for the minimal 5-particle Fock state (a) and for the 7-particle 
Fock states with light quarks q = u, d (b) and with strange quarks (c). The solid curve 
in each case is the charm quark distribution which also serves as the hadron distribution 
for independent fragmentation. The other curves are the probability distributions for 
hadron production by coalescence, including: D~ (dashed), (dot-dashed), (dot- 
dash-dashed) and D~ (dotted). If the shape of the probability distribution is the same 
for any two hadrons (such as the and the A+ in (b)) in a conhguration, it is indicated. 
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Figure 5: Charm hadron production in the intrinsic charm model with a proton pro¬ 
jectile. The probability distributions, {1 /P^){dP^/dx h) , for uncorrelated fragmentation 
and coalescence are given for the minimal 5-particle Fock state (a) and for the 7-particle 
Fock states with light quarks q = u, d (b) and with strange quarks (c). The solid curve 
in each case is the charm quark distribution which also serves as the hadron distribution 
for independent fragmentation. The other curves are the probability distributions for 
hadron production by coalescence, including: D~ (dashed), A+ (dot-dashed), (dot- 
dash-dashed) and D~ (dotted). If the shape of the probability distribution is the same 
for any two hadrons (such as the and the A+ in (b)) in a conhguration, it is indicated. 
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Figure 6: Charm hadron production in the intrinsic charm model with a Ti~ projectile. 
The probability distributions, {1/P^){dP^/dxH), for uncorrelated fragmentation and 
coalescence are given for the minimal 5-particle Fock state (a) and for the 7-particle 
Fock states with light quarks q = u, d (b) and with strange quarks (c). The solid curve 
in each case is the charm quark distribution which also serves as the hadron distribution 
for independent fragmentation. The other curves are the probability distributions for 
hadron production by coalescence, including: D~ (dashed), A+ (dot-dashed), (dot- 
dash-dashed) and D~ (dotted). If the shape of the probability distribution is the same 
for any two hadrons (such as the and the A+ in (b)) in a conhguration, it is indicated. 







Figure 7: Model predictions are compared to the T,~A data of Ref. for (a) D~, 
(b) A+, (c) and (d) The solid and dashed curves represent our full model, with 
the intrinsic charm probability distributions given in eqs. ( |A-1|) - (|A7BD for carbon and 
copper targets respectively. The dot-dashed and dotted curves contrast the results for 
carbon and copper targets respectively with the simplihed model which considers only 
fragmentation from the minimal Fock state and coalescence only from the state with the 
minimum number of partons necessary to produce it. 


29 






















Figure 8: Model predictions are compared to the S A data of Ref. |p for the following 
asymmetries: (a) D~/D'^, (b) D~/Df and (c) A+/A^, as well as our prediction for 
the (d) D~/E^ asymmetry. The solid and dashed curves represent our full model, with 
the intrinsic charm probability distributions given in eqs. ( |A-1|) - (|A-8|) for carbon and 
copper targets respectively. The dot-dashed and dotted curves for carbon and copper 
targets respectively contrast the results with the simplified model which considers only 
fragmentation from the minimal Fock state and coalescence only from the state with the 
minimum number of partons necessary to produce it. 
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Figure 9: Predictions for charm hadron production are given in our full model for S“Cu 
interactions at 650 GeV. The individual xp distributions are given in (a) and (b). All 
cross sections are compared to the leading twist fusion calculation in the solid curve. In 
(a) the hadron distributions are D~ (dashed), A+ (dot-dashed), (dotted) and D~ (dot 
dashed dashed). In (b) the hadron distributions are (dashed), (dot-dashed), 
(dotted) and (dot dashed dashed). Predictions of the asymmetries are given in (c) 
for D~/ (solid) and D~/ (dashed) while the prediction for the asymmetry 

is given in (d). 
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Figure 10: Predictions for charm hadron production are given in our full model for 650 
GeV 7r“Cu interactions. The individual xp distributions are given in (a) and (b). All 
cross sections are compared to the leading twist fusion calculation in the solid curve. 
In (a) the hadron distributions are D~ (dashed), A+ (dot-dashed), (dotted) and Dj 
(dot dashed dashed). In (b) the hadron distributions are (dashed), (dot-dashed), 

(dotted) and (dot dashed dashed). A prediction of the D ~/asymmetry is 
given in (c). 
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Figure 11: Predictions for charm hadron production are given in our full model for pCu 
interactions at 800 GeV. The individual xp distributions are given in (a) and (b). All 
cross sections are compared to the leading twist fusion calculation in the solid curve. In 
(a) the hadron distributions are D~ (dashed), A+ (dot-dashed), (dotted) and D~ (dot 
dashed dashed). In (b) the hadron distributions are (dashed), (dot-dashed), 
(dotted) and (dot dashed dashed). Predictions of the asymmetries are given in (c) 
for D~/ (solid) and D~/ (dashed) while the prediction for the Zl“/A+ asymmetry 
is given in (d). 
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